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ABSTRACT 

We present new SMA and PdBI observations of the intermediate-mass object IRAS 22198+6336 in the 
millimeter continuum and in several molecular line transitions. The millimeter continuum emission reveals a 
strong and compact source with a mass of ~ 5 M© and with properties of Class objects. CO emission shows 
an outflow with a quadrupolar morphology centered on the position of the dust condensation. The CO outflow 
emission seems to come from two distinct outflows, one of them associated with SiO outflow emission. A 
large set of molecular lines has been detected toward a compact dense core clearly coincident with the compact 
millimeter source, and showing a velocity gradient perpendicular to the outflow traced by CO and SiO. The 
chemically rich spectrum and the rotational temperatures derived from CH 3 CN and CH3OH (100-150 K) 
indicate that IRAS 22198+6336 is harbouring one the few intermediate-mass hot cores known at present. 
Subject headings: stars: formation — ISM: individual objects (IRAS 22198+6336) — ISM: lines and bands 
— radio continuum: ISM 



1. INTRODUCTION 

Early stages of massive star formation have numerous 
observational signatures, including molecular masers and 
outflows, sub-millimeter continuum sources, hot molecular 
cores, and hyper/ultracompact HII regions. Hot molecular 
cores (HMCs), first discovered in the vicinity of massive 
protostars (see Kurtz et al. 2000 for a review), are compact 
(<0.1 pc, n > 10 7 cm" 3 ) objects with relatively high tempera- 
tures (7\ > 100 K) that show a very rich chemistry in complex 
organic molecules (CH 3 CN, CH 3 OH, CH3OCHO,. . . ). 

Recently, regions characterized by high temperatures and 
densities have also been detected around low-mass protostars 
(see Ceccarelli 2004 for a review). These so-called hot cori- 
nos share many characteristics with HMCs, although the for- 
mer are ~ 10 2 times less massive. Many questions remain 
unanswered regarding both the high-mass and low-mass hot 
core phase. For example, HMCs are relatively large objects 
(~ 0.1 pc) and harbor clusters of stars, while hot corinos seem 
to be smaller molecular envelopes surrounding a single low- 
mass young stellar object (YSO) or binary. Additionally, al- 
though HMCs are generally understood to be illuminated by 
the newly born high-mass stars embedded within the cores 
themselves (Cesaroni 2005), an alternative scenario has been 
proposed for hot corinos, in which jets and shocks are respon- 
sible for enhancing the temperature and the abundances of 
many species (Chandler et al. 2005). 

Intermediate-mass YSOs (IMYSOs; M* -2-8 M ) are ex- 
cellent targets to study the different hot core heating mech- 
anisms (radiative or shocks), because they share some prop- 
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erties of massive stars, such as their association with clusters 
and their ability to photodissociate and ionize the surround- 
ing gas, and present the advantages of being located closer 
(< 1 kpc) than most of the massive star-forming regions, and 
in regions of less complexity. However, there are only two 
IMYSOs with hot core emission reported in the literature, 
NGC 7129-FIRS 2 and IC 1396 N (Fuente et al. 2005, 2009), 
and few conclusions can be inferred from such a small num- 
ber of intermediate-mass hot cores (IMHCs). 

IRAS 22198+6336 (hereafter 122198) is an IMYSO located 
at a distance of 764 pc (Hirota et al. 2008) with a bolometric 
luminosity of 370 L©. The presence of a strong and com- 
pact submillimeter dust condensation (Jenness et al. 1995) 
with no near-infrared nor mid-infrared emission suggests that 
it is a deeply embedded object, classified as an intermediate- 
mass Class source by Sanchez-Monge et al. (2008). Single- 
dish telescope observations reveal a CO outflow well centered 
on the position of the dust condensation (Zhang et al. 2005). 
Sanchez-Monge et al. (2008) detect slightly resolved partially 
optically thick centimeter emission consistent with a thermal 
ionized wind (or radiojet). Additionally, NH3 and CS dense 
gas emission together with H2O and OH maser emission has 
been detected toward 122198 (Tafalla et al. 1993; Larionov et 
al. 1999; Valdettaro et al. 2002; Edris et al. 2007; Hirota et al. 
2008), all indicative of a very early evolutionary stage. 

In this letter, we present new interferometric continuum and 
molecular line observations that reveal that 122198 is associ- 
ated with an intermediate-mass hot core. 

2. OBSERVATIONS 

We observed the 122198 region with the Submillimeter 
Array (SMA0; Ho et al. 2004) in the 1.3 mm (230 GHz) 
band using the compact array configuration on 2008 June 
11. The phase center was a(J2000) = 22 h 21 m 27 s .617, and 
£(J2000) = +63°51 '42V 18, and the projected baselines ranged 
from 7 kA to 100 LI. The two sidebands of the SMA cov- 
ered the frequency ranges of 219.4-221.4 GHz and 229.4- 
231.4 GHz, with a spectral resolution of ~ 0.5 km s" 1 . Sys- 

7 The SMA is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics, 
and is funded by the Smithsonian Institution and the Academia Sinica. 
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Fig. 1. — IRAS 22198+6336. Left: Blue and red contours: PdBI CO (1-0) intensity map for blueshifted velocities (from -5 to -15 km s _1 with respect to the 
systemic velocity, ^lsr — —11.3 km s _1 ) and redshifted velocities (from +5 to +15 km s _1 with respect to the ^lsr)- Levels start at 5%, increasing in steps of 
15% of the peak temperature, 7200 K km s _1 and 5400 K km s _1 for blue and red respectively. White contours: PdBI SiO (2-1) intensity map for blueshifted 
velocities (from -5 to -15 km s _1 with respect to the ^lsr)- Levels start at 5%, increasing in steps of 30% of the peak temperature, 450 K km s _1 . Synthesized 
beams for CO and SiO images are shown in the bottom right and left corners, respectively. Black contours: SMA 1.3 mm continuum map. Levels are -4, and 4 to 
60 in steps of 8, times 3.92 mJy beam -1 . Grey scale: 4.5 /urn IRAC/Spitzer image. Right: Each of the four panels show in black contours the zero-order moment 
(intensity) map, and in colour scale the first-order moment (velocity) map for four molecular transitions observed with the SMA. Levels are 5%, 45%, and 85% 
the intensity peak 2.67, 3.20, 3.32, 3.97 K km s _1 , for CH 3 CN, HNCO, OCS and CH 3 OH respectively. The (0,0) position corresponds to the SMA phase center. 
The white star marks the position of the centimeter continuum source detected by Sanchez-Monge et al. (2008). 



tern temperatures ranged between 200 and 300 K. The zenith 
opacities at 225 GHz were around 0.15 and 0.20 during the 
track. The FWHM of the primary beam at 1.3 mm was ~ 56". 
Bandpass calibration was performed by observing the quasar 
3C454.3. Amplitude and phase calibrations were achieved 
by monitoring 0019+734 and 1928+739, resulting in an rms 
phase of ~ 30°. The absolute flux density scale was de- 
termined from Callisto and Jupiter with an estimated uncer- 
tainty around 15%. Data were calibrated and imaged with the 
MIRIAD software package. The continuum was constructed 
in the (u,v) domain from the line-free channels. Imaging 
was performed using natural weighting, obtaining a synthe- 
sized beam of 3'.'0 x 2'.'2 with a P.A.=50°, and 1 cr rms of 
80 mJy beam -1 per channel, and 3.9 mJy beam -1 for the con- 
tinuum. 

The IRAM Plateau de Bure Interferometer (PdB0) was 
used to observe the CO (1-0) (115.12 GHz) and SiO (2-1) 
(86.85 GHz) molecular transitions toward 122198. We carried 
out 2-pointing mosaic observations (primary beam^ 56"), in 
2008 and 2009, with the array in the C and D configurations. 
The two spectral setups, tuned at 115.27 GHz (2.7 mm) and 
at 86.85 GHz (3.5 mm), include several molecular transitions, 
which will be presented in a forthcoming paper (Sanchez- 
Monge et al., in prep.). For the CO and SiO we used a cor- 
relator unit of 40 MHz of bandwidth with 512 spectral chan- 
nels, which provides a spectral resolution of ~ 0.3 km s" 1 . A 
number of 320 MHz continuum units were used to image the 
continuum, six at 115 GHz and three at 86 GHz. The typ- 
ical system temperatures for the receivers were ~ 200 K at 
115 GHz and ~ 100 K at 86 GHz. Bandpass calibration was 
performed by observing 3C273 and 3C454.3. Amplitude and 
phase calibrations were achieved by monitoring 0116+731, 
1928+738 and 2037+51 1, resulting in a phase rms of around 



25°. The absolute flux density scale was determined from 
MWC349 and 3C84, with an estimated uncertainty around 
25%. The data were calibrated with the program CLIC, and 
imaged with MAPPING, both parts of the GILDAS0 software 
package. Imaging was performed using natural weighting, 
obtaining a synthesized beam of 370 x 275, with a P.A.=107° 
at 115 GHz, and 475 x 470, with a P.A.=81° at 87 GHz. The 
1 cr rms for the channel and continuum images were 40 and 
0.6 mJy beam" 1 at 115 GHz, and 30 and 0.3 mJy beam" 1 at 
87 GHz, respectivelly. 

3. RESULTS 
3.1. Continuum results 

In Figure Q] (left panel) we show the interferometric mil- 
limeter continuum map toward 122198 (black contours). 
We detected a compact source (with coordinates a = 
22 h 21 m 26 s .78 and 8 = +63°51 , 3775) coincident with the cen- 
timeter source reported by Sanchez-Monge et al. (2008), and 
surrounded by a faint structure extended toward the south- 
west. The source is clearly detected at 3.5, 2.7, and 1.3 mm 
with integrated flux densities of 24 + 6 mJy, 55 + 15 mJy, and 
500 + 75 mJy, respectively, and has a similar morphology at 
the three wavelengths. The deconvolved size of the source is 
370x271, withP.A.=150°, corresponding to 2300 x 1600 AU. 
The millimeter continuum emission is likely tracing the dust 
envelope of 122198. 

The spectral energy distribution (SED) reported by 
Sanchez-Monge et al. (2008) has been improved with the 
new continuum data. We also included the fluxes between 
3.6 yum and 8.0 yum extracted from the Spitzer Space Telescope 
Infrared Array Camera (IRAC) maps (PI: Giovanni Fazio, 
project: 40147). We retrieved the images using the LeoparcQ 



8 The Plateau de Bure Interferometer (PdBI) is operated by the Institut de 
Radioastronomie Millimetrique (IRAM), which is supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). 



9 GILDAS: Grenoble Image and Line Data Analysis System, 
http://www.iram.fr/IRAMFR/GILDAS 

1U http://ssc.spitzer.caltech.edu/warmmission/propkit/spot/ 



The intermediate-mass hot core IRAS 22198+6336 



-10 



-14 



-16 



o IRAS 221 98+6336 (at 760 pc) 
□ NGC7129-FIRS2(at760pc) 




10 10° 10" 

Wavelength (jam) 



10° 



Fig. 2.— Spectral energy distribution of IRAS 22198+6336 (black circles). 
Solid line: sum of the modified blackbody law and the thermal ionized wind. 
Black dashed line: modified blackbody law for the dust envelope with a dust 
emissivity index of fi = 1.7, source radius of 3'/0, dust temperature of 36 K, 
dust mass of 5 M Q , and a dust mass opacity coefficient of 0.9 cm 2 g _1 at 
1.3 mm. Black dotted line: thermal free-free emision from an ionized wind 
with a spectral index of ~ 0.6 at centimeter wavelengths. We estimated the 
total bolometric luminosity, from IR up to mm wavelengths, to be L^ \ =* 
370 L Q for a distance of 0.76 kpc. Red squares: NGC 7129-FIRS 2 data (from 
Eiroa et al. 1998, Crimier et al. 2010, and 2MASS and IRAC catalogues) 
rescaled at the distance of 122198 (760 pc). We estimated Lb i - 480 L©. 
Labels of telescopes refer only to 122198 data. 



software, and evaluated the fluxes using the MOPEXE3 soft- 
ware. The global SED (open circles in Fig. [2]) was fitted by a 
modified blackbody law plus a thermal ionized wind (with a 
spectral index of -0.6 at centimeter wavelengths). Assuming 
a dust mass opacity coefficient given by k v = ko(v/vo)P with 
a:o = 0.9 cm 2 g" 1 at vo = 230 GHz (Ossenkopf & Henning 
1994), the SED is well fitted with a dust emissivity index of 
1.7 ± 0.2, a dust temperature of 36 ± 2 K, and an envelope 
mass of 5 ± 1 M Q . As shown in Fig. El the SED of 122198 
resembles that of Class objects (Andre et al. 1993). 



3.2. Molecular results 

Interferometric 12 CO(1-0) (Fig. ffl left) and 12 CO(2-l) 
(Sanchez-Monge et al., in prep.) maps reveal an outflow with 
a quadrupolar morphology clearly centered on the position of 
the dust condensation. The outflow emission spans a veloc- 
ity range from -25 km s" 1 to +5 km s" 1 for the / = 1 — > 
transition (the systemic velocity is i^lsr — — 11.3 km s" 1 ). The 
quadrupolar morphology of the outflow can be interpreted as 
the walls of a single outflow or as the superposition of two 
bipolar outflows: outflow A in the southwest-northeast direc- 
tion, and outflow B in the northwest- southeast direction. Both 
outflows are centered on the position of the millimeter contin- 
uum source, which could be a millimeter clump harbouring 
at least two sources. In fact, higher angular resolution ob- 
servations (0.4" or 300 AU) with the PdBI at 1.3 mm reveal a 
strong source with indications of a faint extension to the south 
(Palau et al., in prep.). In the following, we will assume the 
two outflow interpretation. We estimated the physical param- 
eters of outflows A and B (see Table [B and found values sim- 
ilar to the outflow parameters found in other IMYSOs (e. g., 
IC 1396 N: Beltran et al. 2002; IRAS 22272+6358 A: Beltran 
et al. 2006). The integrated SiO(2-l) emission, spanning a 
velocity range from -25 to -8 km s" 1 , is shown in white con- 
tours in Fig. [T] This emission coincides with the blueshifted 
lobe of outflow A. No redshifted counterpart is observed for 

11 http://ssc.spitzer.caltech.edu/dataanalysistools/tools/mopex/ 
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Fig. 3. — 230 GHz continuum-free spectra in the image domain from the 
SMA data, for the lower side band (LSB) in the top panel and the upper 
side band (USB) in the bottom panel. The molecular transitions used in the 
rotational diagrams (Fig.|4j are shown in blue and red. 



the SiO emission, however. The detection of SiO in outflow A 
suggests that it is younger than outflow B, because gas-phase 
SiO abundances may decrease with time (Klaassen & Wilson 
2007; Shang et al. 2006). The SiO detection also suggests that 
the CO quadrupolar morphology arises from two distinct out- 
flows, because SiO generally traces material associated with 
the primary jet driving the outflow rather than the walls of 
the outflow cavity (e.g., H211: Palau et al. 2006; Lee et al. 
2007). We also note that the IRAC 4.5 yum band (grey scale 
in Fig. [T] left) shows extended emission ("green fuzzy") as- 
sociated mainly with the blueshifted CO outflow lobes. Such 
4.5 yum extended emission is thought to be associated with 
shocked H2 emission (De Buizer & Vacca 2010). 

The wide-band SMA spectrum toward 122198 (Fig. [3]) re- 
veals a chemically rich dense core with typical hot core fea- 
tures (e. g., Cesaroni et al. 1999; Gibb et al. 2000). Molecular 
line transitions detected with 7b > 0.2 K have been identi- 
fied and imaged. We note that the linewidths observed for 
most of the molecules are 5-8 km s" 1 . Figure [T] (right panel) 
shows the maps of high-density tracers with a compact source 
(2'.'1 x 173, with P.A. = 15°) at the position of the millime- 
ter continuum source. Although the maps of each molecule 
show a source barely resolved, and higher angular resolutions 
should be made, first-order moment (velocity field) maps for 
these high-density tracers (Fig. Q] right) show velocity gra- 
dients in the northwest- southeast direction, perpendicular to 
outflow A, presumable the youngest outflow. Most of these 
molecules trace disks and/or toroids in other massive YSOs 
(e. g., Beltran et al. 2006). However, the presence of outflow 
B perpendicular to outflow A makes the interpretation of these 
velocity gradients uncertain. Assuming that these gradients 
trace rotation, the dynamical mass has been calculated from 
the expression Md yn = i% ot RvotIG sin 2 /, where f; rot is the ve- 
locity estimated from the gradient, R m the radius of the core, 
G the gravitational constant, and i the inclination angle as- 
sumed to be 90° (edge-on). In our case, with v YOt ^ 2 km s" 1 
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Fig. 4. — CH3CN (blue squares) and CH3OH (red circles) rotational dia- 
grams. Open symbols and dashed lines refer to the optically thin approxi- 
mation, while filled symbols and solid lines refer to the opacity which yields 
the best linear fit. For the optically thick fit we give the values in the figure, 
and adopt them in the main text. For the optically thin fit, T wt = 160 ± 30 K 
and AW = (1-7 ± 0.3) x 10 14 cm" 2 for CH3CN, and T mt = 220 ± 60 K and 
Af mol = (4.6 + 1.6) x 10 16 cm -2 for CH3OH. Vertical left and right axes refer 
to CH3CN and CH3OH, respectively. 

and R m ^ 1 " for the four molecules shown in Fig. Q] right, 
the dynamical mass is ~ 3.5 M©. Additionally, following the 
models of Palla & Stahler (1993), and assuming a luminosity 
for the object of ~ 370 L© (Fig. [2]), the mass of the embedded 
YSO is about 4-5 M Q . Thus, the mass of the internal YSO 
estimated from both methods is ~ 4 M Q , which is similar to 
the dust envelope mass (~ 5 M©), as is characterise of Class 
objects. 

Because of the wide SMA band, one can observe several 
CH3CN and CH3OH transitions simultaneously. Following 
the rotational diagram method (Goldsmith & Langer 1999; 
Araya et al. 2005), which assumes that all molecular levels 
are populated according to the same excitation temperature, 
we can derive this temperature, T ToU and the total column den- 
sity, N mo \, for these molecules. We constructed the rotational 
diagrams following two methods. First, we assumed optically 
thin emission (open symbols in Fig. [4]). Second, we estimated 
the opacity for each transition (following Goicoechea et al. 
2006; see also Girart et al. 2002) and we adopted the value 
that yielded the best linear fit (filled symbols in Fig. [4]). The 
resulting opacities for the different CH3CN and CH3OH tran- 
sitions range from 0.3 to 3.5 and from 0.2 to 6.5, respectively. 
In the discussion we will use the results of the latter fit, which 
are T mt = 100 ± 30 K and N mo \ = (3.6 ± 0.4) x 10 14 cm" 2 
for CH 3 CN, and T m = 150 ± 40 K and N mol = (8.6 ± 1.3) x 
10 16 cm" 2 for CH3OH. The high temperatures, together with 
the chemically rich spectrum, are clear evidence that 122198 
is an intermediate-mass hot core. 

4. DISCUSSION AND CONCLUSIONS 

Our interferometric continuum and molecular line observa- 
tions toward 122198 reveal a compact and strong millimeter 
source driving a quadrupolar outflow, which is associated with 
a compact dense core with hot core signatures. At infrared 
wavelengths, IRAC images reveal a single source coincident 
with the hot core and with "green fuzzies" at 4.5 //m. The 
high rotational temperature (100-150 K) estimated from the 
CH3CN and CH3OH rotational diagrams, and the chemically 
rich spectrum, make 122198 one of the few IMHCs known at 
present. There are only two other IMHCs reported in the liter- 
ature: NGC 7129-FIRS 2 (Fuente et al. 2005), and IC 1396 N 
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a References are: 122198: Sanchez-Monge et al. 2008; this work (the 
two values in outflow parameters correspond to outflows A and B, 
respectively); NGC 7129-FIRS 2: Fuente et al. 2001, 2005, Crimier et 
al. 2010 (centimeter data from the VLA project AR304, and outflow 
parameteres: Sanchez-Monge, priv. com.); IC 1396 N: Beltran et al. 
2002; Neri et al. 2007; Fuente et al. 2009. Upper limits for IC 1396 N 
are due to the cluster properties (see text). 

b Estimated following Fuente et al. (2005), assuming T = 100 K and a 
size of 271 x 173, which is the deconvolved size of the CH3CN core. 

(Neri et al. 2007; Fuente et al. 2009). Table Q] lists the main 
properties of the three IMHCs known to date. 

IC 1396 N is associated with a single IRAS source, 
IRAS 21319+5802. However, high-angular resolution mil- 
limeter images reveal a cluster of sources, only one of which 
is associated with the IMHC (Neri et al. 2007). Moreover, 
while Spitzer/IRAC observations reveal a counterpart for mil- 
limeter source 41.73+12.8, there is no counterpart for the 
IMHC millimeter source (41.86+1 1.9; Choudhury et al. 2010; 
Neri et al. 2007). All this makes it difficult to assess the IMHC 
contribution to the IRAS and IRAC fluxes, and hence we re- 
frain from comparing IC 1396 N with 122198. 

The other IMHC, NGC 7129-FIRS 2, has a bolometric lu- 
minosity similar to 122198, has a single associated source 
dominating the millimeter emission (Fuente et al. 2005) and 
also is clearly associated with a source detected in IRAC 
and MIPS bands (Sanchez-Monge, priv. com.). Eiroa et al. 
(1998) study the far-infrared emission, and find one HIRES 
source, which is coincident with the MIPS source (Crimier et 
al. 2010). The association of the millimeter, far-infrared and 
mid-infrared source with the IMHC is clear, as in the case of 
122198 and, similarly, "green fuzzies" in the 4.5 fim IRAC 
image appear associated with the molecular outflow, which 
is also quadrupolar (Fuente et al. 2001) and has physical pa- 
rameters similar to the outflows in 122198 (Table [T]). Regard- 
ing the centimeter emission, the upper limit for NGC 7129- 
FIRS 2 is consistent with the centimeter emission found to- 
ward 122198 (Table [T]). Furthermore, Figure [2] shows that 
both IMHCs have quite similar Class 0-like SEDs. Thus, 
NGC 7129-FIRS 2 and 122198 have very similar properties 
concerning clustering, centimeter and outflow emission, and 
have a similar SED. Moreover, the hot core temperature, T voU 
and CH3CN abundance, Xch3Cn, are similar for these two 
IMHCs. At this point, our data do not allow us to distin- 
guish between the proposed hot core scenarios (radiative or 
shocked), because the two IMHCs have similar luminosities 
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and outflow/centimeter properties. A detailed study of dif- 
ferent molecular abundances (enhanced by the radiation field 
or by shocks) will be presented in a forthcoming paper, and 
should clarify the formation scenario of the hot core. 

In summary, the properties of the newly discovered IMHC, 
122198, and its comparison with the other known IMHCs help 
us better establish the characteristics of the hot core phase 
for IMYSOs. We have found that the IMHC phase is co-eval 
with mass ejection (as shown by outflow, "green fuzzies", and 
faint centimeter emission), presents Class 0-like SEDs, has 
hot core temperatures ~ 100 K, and Xch 3 cn - 10~ 9 - A de- 
tailed comparison of 122198 and NGC7129-FIRS 2 and the 
discovery of new IMHCs will provide fundamental tools to 
understand the hot core phase of star formation. 



We thank the anonymous referee for his/her useful comments. 
A. S-M. and A. P. are grateful to Esteban Araya for his help in ana- 
lyzing the data, to Francesco Fontani and Asuncion Fuente for kindly 
providing complementary data, and also to Rob Gutermuth for his 
help with Spitzer data. A. S-M., A. P. and R. E. are supported by 
the Spanish MICINN grant AYA2008-06189-C03 (co-funded with 
FEDER funds). A. P. is supported by a JAE-Doc CSIC fellowship 
co-funded with the European Social Fund. This research has been 
partially funded by Spanish MICINN under the ESP2007-65475- 
C02-02 and Consolider-CSD2006-00070 grants. S.K. is partially 
supported by DGAPA-UNAM grant IN101310. The National Radio 
Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Univer- 
sities, Inc. 



REFERENCES 



Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 
Araya, E., Hofner, P., Kurtz, S., Bronfman, L., & DeDeo, S. 2005, ApJS, 157, 

279 
Beltran, M. T, Girart, J. M., & Estalella, R. 2006, A&A, 457, 865 
Beltran, M. T, Girart, J. M., Estalella, R., Ho, P. T. P., & Palau, A. 2002, ApJ, 

573, 246 
Ceccarelli, C. 2004, Star Formation in the Interstellar Medium: In Honor of 

David Hollenbach, 323, 195 
Cesaroni, R. 2005, Massive Star Birth: A Crossroads of Astrophysics, 227, 

59 
Cesaroni, R., Felli, M., Jenness, T, Neri, R., Olmi, L., Robberto, M., Testi, 

L., & Walmsley, C. M. 1999, A&A, 345, 949 
Chandler, C. J., Brogan, C. L., Shirley, Y. L., & Loinard, L. 2005, ApJ, 632, 

371 
Choudhury, R., Mookerjea, B., & Bhatt, H. C. 2010, A pJ, 717, 1067 
Crimier, N., Ceccarelli, C, Alonso-Albi, T, et al. 2010, arXiv: 1005.0947 
De Buizer, J. M., & Vacca, W. D. 2010, AJ, 140, 196 
Edris, K. A., Fuller, G. A., & Cohen, R. J. 2007, A&A, 465, 865 
Eiroa, C, Palacios, J., & Casali, M. M. 1998, A&A, 335, 243 
Fuente, A., Neri, R., Martin-Pintado, J., Bachiller, R., Rodriguez-Franco, A., 

& Palla, F. 2001, A&A, 366, 873 
Fuente, A., Neri, R., & Caselli, P. 2005, A&A, 444, 481 
Fuente, A., Castro-Carrizo, A., Alonso-Albi, T, et al. 2009, A&A, 507, 1475 
Gibb, E., Nummelin, A., Irvine, W. M., Whittet, D. C. B., & Bergman, P. 

2000, ApJ, 545, 309 
Girart, J. M., Viti, S., Williams, D. A., Estalella, R., & Ho, P. T. P. 2002, 

A&A, 388, 1004 



Goicoechea, J. R., Pety, J., Gerin, M., Teyssier, D., Roueff, E., Hily-Blant, P., 

& Baek, S. 2006, A&A, 456, 565 
Goldsmith, P. E, & Langer, W. D. 1999, ApJ, 517, 209 
Hirota, T, Ando, K., Bushimata, T, et al. 2008, PASJ, 60, 961 
Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 
Jenness, T, Scott, P. E, & Padman, R. 1995, MNRAS, 276, 1024 
Klaassen, P. D., & Wilson, C. D. 2007, ApJ, 663, 1092 
Kurtz, S., Cesaroni, R., Churchwell, E., et al. 2000, Protostars and Planets IV, 

299 
Larionov, G. M., Val'tts, I. E., Winnberg, A., Johansson, L. E. B., Booth, 

R. S., & Golubev, V. V. 1999, A&AS, 139, 257 
Lee, C.-E, Ho, P. T. P., Palau, A., et al. 2007, ApJ, 670, 1 188 
Neri, R., Fuente, A., Ceccarelli, C. et al. 2007, A&A, 468, L33 
Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 
Palau, A., Ho, P. T. P., Zhang, Q., et al. 2006, ApJ, 636, L137 
Palla, E, & Stahler, S. W. 1993, ApJ, 418, 414 
Sanchez-Monge, A., Palau, A., Estalella, R., Beltran, M. T, & Girart, J. M. 

2008, A&A, 485, 497 
Shang, H., Allen, A., Li, Z.-Y, Liu, C.-E, Chou, M.-Y, & Anderson, J. 2006, 

ApJ, 649, 845 
Tafalla, M., Bachiller, R., & Martin-Pintado, J. 1993, ApJ, 403, 175 
Valdettaro, R., Palla, E, Brand, J., Cesaroni, R., Comoretto, G., Felli, M., & 

Palagi, F. 2002, A&A, 383, 244 
Zhang, Q., Hunter, T. R., Brand, J., Sridharan, T. K., Cesaroni, R., Molinari, 

S., Wang, J., & Kramer, M. 2005, ApJ, 625, 864 
Zhang, Q., Hunter, T. R., Beuther, H., Sridharan, T. K., Liu, S.-Y, Su, Y.-N., 

Chen, H.-R., & Chen, Y 2007, ApJ, 658, 1152 



